Nanowire waveguide made from extremely anisotropic metamaterials 



00 

o 
o 

<N 

Oh: 

< 

OO 



C/3 

O 

-(— > 
Oh" 

q 
o 

Oh. 



> 

in 

o 

00 

o 



X 



Y. J. Huang, W. T. Lu,* and S. Sridhar 

Department of Physics and Electronic Materials Research Institute, 

Northeastern University, Boston, Massachusetts 02115, USA 

(Dated: April 28, 2008) 

Exact solutions are obtained for all the modes of wave propagation along an anisotropic cylindrical 
waveguide. Closed-form expressions for the energy flow on the waveguide are also derived. For 
extremely anisotropic waveguide where the transverse permittivity is negative [e^ < 0) while the 
longitudinal permittivity is positive {e\\ > 0), only transverse magnetic (TM) and hybrid modes will 
propagate on the waveguide. At any given frequency the waveguide supports an infinite number 
of eigenmodes. Among the TM modes, at most only one mode is forward wave. The rest of them 
are backward waves which can have very large effective index. At a critical radius, the waveguide 
supports degenerate forward- and backward- wave modes with zero group velocity. These waveguides 
can be used as phase shifters and filters, and as optical buffers to slow down and trap light. 

PACS numbers: 



I. INTRODUCTION 

Since the realization of negative refraction [1] in mi- 
crowaves [2], there is renewed and intense interest in 
electromagnetic metamaterials. Negative refraction has 
added a new arena to physics, leading to new concepts 
such as perfect lens [3, 4], superlens [3, 5, 6], and focusing 
by plano-concave lens [7, 8] . Negative refraction has sub- 
sequently been realized in microwaves [9-13], THz waves, 
and optical wavelengths [14-16], in metamaterials made 
of wire and split-ring resonators [17] or photonic crystals 
[18-20]. 

Metamaterials are artifically fabricated structures pos- 
sessing certain desirable properties which are not avail- 
able in natural materials. Metamaterials can have double 
negative index [21] or single negative index. Metamateri- 
als can be periodic, such as photonic crystals [22]. They 
can also be non-periodic, such as the materials for cloak- 
ing [23]. They can also be made to be anisotropic and 
have indefinite index [24-26]. Indefinite index matemate- 
rials can be used to make hyperlens [27, 28]. This range 
of properties opens an infinite possibilities to use mate- 
materials in frequencies from microwave all the way up 
to the visible. 

Wave propagation in waveguide of nanometer size [29] 
has unique properties. In this paper, we consider wave 
propagation along anisotropic nanowires. In the case 
where the transverse permittivity is negative while the 
longitudinal one is positive {e_\_ < 0, ^n > 0), these in- 
definite index waveguides can support both forward and 
backward waves. High effective index can be obtained 
for these modes. These waveguides can also support de- 
generate modes which can be used to slow down and trap 
light. 

In Sec. II, we derive the formulas for all the modes 
on the anisotropic cylinders. Exact solutions for all the 



modes and closed-form expressions for the energy flow 
will be obtained. Possible zero net-energy flow modes 
will also be discussed. The situation for trapping light 
is presented in Sec. III. In Sec. IV, we propose the re- 
alization of nanowires made of indefinite index medium, 
which is confirmed in finite- difference time-domain simu- 
lations. We conclude in Sec. V with possible applications 
for these anisotropic nanowires. 



II. WAVE PROPAGATION AND ENERGY 

FLOW ON ANISOTROPIC CYLINDRICAL 

WAVEGUIDES 

We consider wave propagation on a cylindrical waveg- 
uide. The axis of the waveguide is along the z-direction 
as shown in Fig. 1. The waveguide is nonmagnetic and 
has an anisotropic optical property 



£y — £t 7^ £z' 



(1) 




s, >0 



FIG. 1: (Color online) An anisotropic cylindrical waveguide 
with axis along the 2;- axis. The longitudinal permittivity is 
positive while the transverse permittivity is negative. 

The waves will propagation along the cylinder axis 
with 
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E = Eoe'(''^-'"*) , H = Hoe^C'^-"*) . 



(2) 



Here P is the propagation wave number along the waveg- 
uide. 

Due to the symmetry of the waveguide, ah the field 
components can be expressed in terms of the longitudinal 
components Ez and Hz. In the polar coordinate system, 
one has for the fields inside the waveguide r < a with a 
the radius. 



Er = i 



^-—(^fSdrEz^ko^d^Hz), 



stk^ - (3- 



Hr = i 



etk§ - p'^y r 

^4> = «^3 — -r^(stkodrE,+l3-d^H,). (3) 

Here ko is the wave number in the vacuum. 

The wave equations for the longitudinal components 
inside the waveguide are 

(d^ + dl)E, + £,(fc2 - l3^let)E, = 0, 
{dl + dl)H, + {stkl - /?2)J?, = 0. (4) 

The waveguide is free-standing in air, so the wave equa- 
tions for r > a are given by the above equations with 
the permittivity replaced by unity. The solutions are ex- 
pressed in terms of the Bessel functions of various kinds 



Ez = AJn(i^r)e^"^, r < a, 
= Ci^n(A^or)e^^^, r>a 



and 



Hz = BIn{^r)e 



in(p 



r < a, 



= DK„{Kor)e''"l' , r > a. 



(5) 



(6) 



The coefficients will be determined by matching the 
boundary conditions. Here 



K = ^z^kl-P^/st 



^ = V^^ -^tkl, 



kq = y^2 _^2^ 



(7) 



We only consider the extremely anisotropic case such that 
the longitudinal permittivity is positive while the trans- 
verse permittivity is negative 



Sz > 0, et < 0. 



(8) 



One can see that due to the anisotropic nature of the 
waveguide, Hz and Ez inside the waveguide will have 
completely different behaviors. 

The continuity of Ez and Hz at the interface r = a 
gives 



C Jn{Ka) D In{i^a) 



A Kn{i^oct)' B Kni^oci)' 



(9) 



The continuity of E^jy at the interface gives 



injS- 



kQa^nQK? A Jn{Ka) 
with the following defined functions 

K{X) ^ Kn-l{x) 
xKn 

In-l{x) 



gn{i^oci) ^hn{t^a) (10) 



9n{x) 



xKn{x) xKn{x) x'^ ' 



hn{x) 



Ini^) 



n 



Xln{x) Xln{x) x' 

The continuity of i^^ at the interface gives 



(11) 



in(3 



with 



Kj Hvn 



AJn{Ka) 



kQa^n^K? B In{t^o) 



gnii^oa) - Szfn{Ka) (12) 



fn{x) 



J'^{x) Jn-l{x) 



n 



(13) 
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Thus we obtain the equation for all the modes 

= n^ {K,oa)~'^ — {K,a)~'^ {i<iQa)~'^ — et{i^a)~'^ (14) 

For wave propagation on the cylindrical waveguide, the 
components of the Poynting vector are 



Of 

Sch 



1 

47r 



{ErH^ - Eci^H*), 



47r 



{EzH;-e^h:), 



47r 



{EzH; - ErHl). 



(15) 



The physical Poynting vector is given by 3?S. 

The total energy fiow along the waveguide is the sum 
of energy fiow inside and outside the waveguide 



P — P 

^z — J^^ 



m I Dout 



(16) 



with 



-I 



27r / Szvdr, P°^* = 27r / Szvdr. (17) 



/ 

da 



Following Ref. [30], the total energy fiow is normalized 
as 



(n) 



-)\w \_ pout 



I pout I 



(18) 



Thus one has -1 < (P^) < 1. 

In the following, we discuss different modes in detail. 



TE modes 



For the transverse electric (TE) modes, Ez = 0. The 
longitudinal magnetic field is given by Eq. (6). One has 

Ed) = i^rdrHz = i — Bli(hir), r < a, 

= i^drHz = -i—DKiUr), r > a. (19) 

The continuity of £^0 at the interface requires that 

ho{K.a) + go{K.oa) = 0. (20) 

For materials without loss, each term on the left side is 
positive, thus there is no solution. The waveguide does 
not support TE modes. This is exactly like that of a 
metallic wire which does not support TE surface waves 
since current must fiow along the waveguide. 

Only when e^ > 1, the waveguide will support TE 
modes, like an ordinary dielectric fiber. 



B. TM modes 

For the transverse magnetic (TM) modes, Hz = 0. The 
longitudinal electric field is given by Eq. (5). One has 

H^ = -iez-^AJi{Kr), r < a, 

= i — CKi{f^or), r>a. (21) 

The continuity of H^f) leads to the equation 

ezfo{Ka)=go{^oa). (22) 

The solutions to this equation give all the TM modes. 

1. Band structure of TM modes 

We first consider the solutions for fixed and real values 
of Sz and St. This is normally associated with a fixed /cq. 
It is convenient to consider solution in the form of Ka or 
the reduced radius koa as a function of tzoa. The wave 
number along the waveguide can be obtained through 
/3 = {k^ -\- /^o)"^^^- Before we seek general solutions, it 
is better to consider the solutions in certain limits to 
reveal some important features of the TM modes on the 
anisotropic waveguide. 

For the TM modes close to the light line, /t:o<^ ^ 0, one 
has 



Ka 



^0,1 



XO,' 



-{i^oaf{\n— +7). 



Here we have used Ko{x) = — ln(x/2) — 7 for small au- 
gument with 7 the Euler constant. For complex St with 
SRst < and ^et > 0, the real and imaginary parts of /? of 



the allowed modes will have the same signs. These modes 
are forward waves, similar to that of an ordinary optical 
fiber. We note that close to the light line, the property 
of the TM modes of the anisotropic waveguide is similar 
to that of an isotropic fiber with s = 1 -\- Sz{l — s^^). 

In the limit of long wavelength or small waveguide ra- 
dius, koa <C 1, Eq. (22) is reduced to 

£zfo{i^oa/r]) = go{i^oa) (23) 



with T] = y—etjsz. This equation gives an infinite num- 
ber of solutions K^a = ^o,m with m = 1, 2, 3, • • •. This in- 
dicates that the anisotropic waveguide supports infinite 
number of propagating modes, no matter how thin the 
waveguide is. For n^a -^ 00, since go{i^oCi) — {tzoCi)~^ -^ 
0, one has ^rn — VXi,m- Here Xn,m is the m-th zero of 
Jn{x) away from the origin. For the Tn-th TM band, 
one has < tzoa < ^o,m- The m-th band starts with 
koa = xo^ral \/^z — ^z/^t whcu K^od = and ends at 
koa = when noa = ^o,m- The modes with no ^ ko 
have dikoa) / d{noa) < and are backward wave. It will 
be obvious if we include small imaginary part in St with 
^St > 0. The equation will give /3 with the real and 
imaginary parts having opposite signs. The energy fiow 
is opposite to the phase velocity, which will be discussed 
later. 

For arbitrary values of /t:o«, the solution must be sought 
numerically. Since the right-hand side of Eq. (22) is 
always positive, the solution requires that Ji{Ka) and 
Jo{Ka) have different signs. For the m-th band, since 
< Koa < ^o,m with ^o,m the solutions of Eq. (23), 
one has xo,m ^ Ka < ^o,m/v < ^i,m- For each tzoCi 
value, the Ka value can be searched within [xo,m5^i,m] 
to satisfy Eq. (22). Once the corresponding iCa is found, 
the reduced radius can be obtained as 



koa : 



V{-et/ez){Kay-{Koay 



For the Tn-th band, the corresponding transverse electric 
field Ez will have m nodes. The band structure for a 
waveguide with et = —3 and e^ = 2 is shown in Fig. 2. 
The effective index of the waveguide neff = /3/ko is also 
evaluated and plotted in Fig. 3. 

Unlike an ordinary fiber where for each band, 
d{koa) / d{noa) > 0, each TM band of the anisotropic 
waveguide starts with d{koa) / d{Koa) > for small nod 
or near the light line. At certain value of nod or koa 
which is marked in Fig. 2, d{koa) / d{noa) — 0. Further 
increasing Koa results in d{koa) / d{Koa) < 0. The band 
ends at a finite Koa = ^o,m- Immediately below the point 
koa where d{koa) / d{Koa) = 0, each band has two modes 
with opposite signs of d{koa) / d{Koa) . One mode is for- 
ward wave and the other backward wave. This will be 
discussed later in the paper. 

We next consider a waveguide of a fixed radius a with 
the following permittivity 



et 



\{l+e. 



S/^o)' 



pi 

2\ I\u2( 



= 2ea{H-k'^)/[kt{l+ea) 



(24) 




FIG. 2: Band structure of the guided TM modes on an 
anisotropic waveguide of radius a with st = —3 and e^^ = 2. 
Open circles denote the degenerate points of forward-wave 
and backward- wave modes. The dashed hues are for a dielec- 
tric waveguide with e = 1 -\- €z{'\- — £^'^)- 




Reduced radius k a 



FIG. 3: Effective index rieff as a function of the reduced radius 
koa for the first three bands of TM modes on an anisotropic 
waveguide with £t = —3 and e^ = 2. Open circles denote the 
degenerate points of forward- wave and backward- wave modes. 



Here Sa and kp are positive constants. The realization 
of this property will be discussed later in Sec. IV. If 
ko < kp/{l + Sa)"^^^, one has et < and Sz > 0. The 
band structure of the TM modes on this waveguide is 
obtained by numeric means and plotted in Fig. 4 with 
the corresponding effective index in Fig. 5. For this 
waveguide, there is no cutoff of hzoa for each band. This 
is because that as /cq ^ 0, e^ — '^^a and et -^ — oo, thus 
V = \/ —^tl^z -^ oo. The cutoff ^o,i 



: r]Xx: 



oo. 



We point out that the Fade approximant for the func- 
tion /o(^) = —x~^ J\{x) I J{){x) can be used to obtain 
good estimate of the solutions. This will be discussed in 
the Appendix. 

If e^ > and Sz <^^ the waveguide also supports TM 
modes. The details will not be presented here. 




FIG. 4: Band structure of the guided TM modes on an 
anisotropic waveguide of radius a with St and Sz given by 
Eq. (24) with a = lO/Zcp and Sa = 2.25. 
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FIG. 5: Effective index rieff as a function of the reduced wave 
number /coa for the first three bands of TM modes on an 
anisotropic waveguide with St and Sz given by Eq. (24) with 
a — lO/Zcp and Sa — 2.25. 



2. Energy flow of TM modes 

Within the waveguide, one has Ez = AJo{Kr)^ E^ = 
i{ezP/etK)AJi{Kr), and H^ = i{ezko/K)AJi{Kr). So 
the Foynting vector component along the axis of the 
waveguide is 



p 



etkoo? 



Ji{Kr) 



JiiKa) 



(25) 



Here we set the coefficient A = K/[ezkoaJi{Ka)]. Since 
Et < 0^ the energy flow inside the nanowire is always 
opposite to the phase velocity. 

For the field in the air r > a, one has Ez = Ci^o(^o^), 
Er = i{Pf^o)CKi{f^or), and H^ = i{ko/f^o)CKi{f^or), 
thus 



P 



Ki{Kor) 



Here the coefficient C 



k^o? Ki{tvoa) 

i^o/[koaKi{i^oa)]. 



(26) 



For the TM modes, one has 



3. Forward-wave and backward-wave TM modes 



pm 

pout 



/? 



P 
2koa? 



MKr) 



Ji{Ka) 
KiiKor) 2 



rdr, 



Ki{Koa) 



rdr. 



(27) 



The above integrals can be carried out and more compact 
expressions for the energy flow can be obtained as 



pm _ 



f3 



Aetko{Ka)nf§{Ka) fo{Ka) 
4koeif§{Ka) StKa ' 



P 



{Kaf 



4fco(Koa)2 Vgl{Koa) go{Koa) 



{Koaf 



P 



JoiKoa) 



4fco(/o('*offl) Koa 



(28) 



Here ^o(^) ^^^ /o(^) ^^^ ^^^ derivatives of ^o(^) and 
/o(x), respectively. 

For convenience, we set ^^ > throughout the paper. 
Since ^o(^) < ^ and f^ix) < 0, one has P^" < and 
pout > 0. In this convention, if (P^) > 0, this indicates 
that the energy flow and the phase propagation are in the 
same directions and the mode is a forward-wave mode. 
Otherwise (P^) < 0, the group velocity and the phase 
velocity are in the opposite direction and the mode is a 
backward-wave mode. The normalized energy flow for 
TM modes on a waveguide with et = —3 and e^ = 2 is 
shown in Fig. 6. We note that for some portion of the 
bands the value of (P^) is negative and thus these modes 
are backward waves. 




Reduced radius k^a 

FIG. 6: (Color online) Normalized energy flow {Pz) for the 
first four bands of the TM modes on an anisotropic waveguide 
with £t = — 3 and e^ = 2. Here we set /3 > 0. Open circles 
denote the degenerate points where (Pz) = 0. 



There are three ways to determine whether a mode 
is a forward wave or backward wave. One is through 
the sign of the derivative d{kQa) / d{i<ioa) . From the band 
structure, we have already notice that for the modes 
near the light line, d{koa) / d{Hioa) > 0. These modes 
are forward waves. For large hzoa or small /coa, one has 
d{koa) / d{f<ioa) < 0, these modes are backward waves. 
From the band structure shown in Fig. 2, the majority 
of the modes are backward waves. 

The second way is through the sign of (Pz)- For the 
TM modes when Koa -^ oo, one has fo{Ka) -^ since 
go{hZoa) -^ 0. One thus has Ka -^ xo^rn- This so- 
lution leads to the divergence of P^" which is negative 
and the vanishing of P°^* which is positive, subsequently 
(Pz) ^ — 1, these modes are all backward waves. Corre- 
spondingly, one has d{koa) / d{hZoa) -^ — oo for hzoa -^ oo. 
This is evident from the band structure shown in Fig. 2. 

In the following, we prove that d{koa) / d{f<ioa) > leads 
to (Pz) > and vice versa. We consider the derivative 



d{Ka) _ 
d{nQa) 

Thus one has 

d{Ka) 



SfKa 



/^o« + (1 - St)koa 



d{koa) 



d^KQo) 



Szt^Qd 



d{Koa) StKa 



>0 if 



d{koa) 
d{K,oa) 



>0. 



On the other hand, according to the eigen equation 
£zfo{Ka) = go{K.oa), one has 

d{Ka) g'oii^oci) 



d{Koa) SzfoiKa)' 

Making use of the above expression, one arrives at the 
inequality 



elf^{Ka) _ g'^^t^Qo) 
£tKa n^a 



>0, 



and subsequently 

pm ^ pout ^ Q^ 

Similarly one has 

pin ^ pout < Q 



if 



if 



d{koa) 
d{nQa) 



d{koa) 
d{nQa) 



>0. 



<0. 



(29) 



(30) 



So the condition for P^ = can be allocated from 
the band structures as shown in Fig. 2, 4 when 
d{koa) / d{f<ioa) = 0. 

The third way to determine whether a mode is a for- 
ward or backward wave is through the relative sign of 
the real and imaginary parts of /3 if dissipation is in- 
cluded. For example we consider et = —3 + 0.05i and 
Sz = 2. At koa = 1.6, the wave numbers of the first three 



eigenmodes are f3a = ±(1.7112 + 0.0067i), ±(2.7250 - 
0.0397i), ±(7.5756 — 0.0676i). Since the free space wave 
length is A = 3.927a, this is a subwavelength waveguide. 
For the TM modes, exept for the first mode, ah the 
other modes are backward-wave modes since for those 
modes 3?/? and ^P have different signs. The normal- 
ized energy flow is (P^) = 0.5151 - 0.0020^,-0.4002 - 
0.0058Z, -0.8760-0.0078 for the above three modes, re- 
spectively. Here we set 3?/3 > 0. The field and Poynting 
vector profiles are plotted in Fig. 7. 



C. Hybrid modes 

The modes with both E^ y^ and H^ y^ are called 
hybrid modes. Their dispersions are contained in the 
solutions of Eq. (14) with n 7^ 0. We recast the equation 
in the following form 

e.UKa) = UO g^iO+h^iy) ' ^''^ 

Here we use the notati on ^ = /^pg and y = na. Note that 
Ka = y/f] with r] = y^-St/e^. 
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FIG. 7: (Color online) The longitudinal electric field and 
Poynting vector for the first three TM mode on an anisotropic 
waveguide of radius a with St = — 3±0.05i and s^^ = 2 at koa = 
1.6 with (a) f3a = 1.7112 + 0.0067i, (b) pa = 2.7250 -0.0397i, 
(c) pa — 7.5756 — 0.0676i. The imaginary parts ^Ez and ^Sz 
are small and not plotted. 



There is an interesting feature of the modes on the 
a nisotropic wave guide. At a fixed /cq, for a < am = 
y/—et/sz{l — St)xo^rn/ko^ the 7Ti-th band TM modes are 
backward waves. If the radius a > am^ the waveguide 
supports two TM modes for the m-th band, one forward 
and one backward. At a = ac^ these two modes be- 
come degenerate and the total energy flow is zero. This 
can be seen in Fig. 2, 3, 6 where degenerate points are 
marked. Further increasing the radius, the waveguide 
will no longer support the 7Ti-th band. The critical ra- 
dius ttc is located such that (P^) = 0, d{koa) / d{i<ioa) = 
or duef^/da = 00. These degenerate modes can be used 
to slow down and even trap light. This will be discussed 
in the next section. 



1 . Band structure of hybrid modes 

At a fixed wavelength A or wave number /cq, £t, ^z are 
constant. Since Ka = y/77, if we use ^ = n^a as a free 
parameter, the eigen equatio n gives a single valu e of Ka 
or y for each ^. Since y = ^J^^ ± (1 — ^t)(^o<^)^, so the 
eigen equation actually gives the reduced radius k^a for 
each ^ = K^a. 

In the limit of long wavelength or small waveguide ra- 
dius, k^a -^ 0, Eq. (31) is reduced to 



£zfn{^/v) =9n{0' 



(32) 



This will give a discrete set of solutions hcoa = £,n,m foi" 
each n. For 77 and Sz not very small, the values £,n,m 
can be obtained approximately by making use of the 
asymptotic behavior gn{x) ~ x~^ for x -^ 00 and the 
Fade approximant of fn{x)^ which will be discussed in 
the Appendix. The anisotropic waveguide supports in- 
finite number of hybrid modes, no matter how thin the 
waveguide is. 

Close to the light line, ^ ^ 0, the eigen equation can be 
simplified. We consider the hybrid modes with different 
n separately. 

For n = 1, one has gi{x) 0^ x~^ — ln(x/2)— 7 for x ^ 0, 
thus Eq. (31) becomes 



Ji{Ka) 



-2(ln^ 
^ 2 



•7 



hi{y) 



9 

y 



(33) 



with ^ <C y. One has the solution Ka -^ xi^rn for ^ ^ 0. 
Note that throughout this paper, we use Xn,m to de- 
note the m-th zero of Jn{x). Also that Xn,o = for 
n > 1. Since xi^o = 0, special care should be taken 
for solutions Ka ~ 0. Making use of the asymptotic 
behaviors hn{x) o:^ nx~'^ ± l/[2(n ± 1)] and /n(^) — 
nx~'^ — l/[2(n ± 1)] for X ^ with n > 1, one gets 



Ka 



1 



l^st 



7?Vln(e/2)±7-(^.±l)/8 



(34) 



From this expression, one can see that only if e^ < —1, 
there will be solution for Ka ~ when ^ ^ 0. The above 
expression gives the dispersion of the first band with n = 
1. The m-th band will start with Ka = Xi m-i and end 



with Ka = £,i^m/V' One has ^i^i < 77x1^1. However if 
— 1 < £t < O7 there wih be no solution for Ka < xi^i. 
All the allowed modes will have finite Ka. For the ?77.-th 
band, one has xi^rn ^ Ka < ^i,m/^ with ^1^1 > 7^X1^1. 
The first band starts with Ka = xi^i. For any xi^rn > 0, 
one has 



Ka :^ xi^r, 



StXl^n 



(^^i,m)'[-2(ln(e/2) +7) + hi{r]xi,m)] ^ St - 

(35) 
The hybrid modes near the light line with n = 1 are all 
forward waves. 

For n > 3, one has the following asymptotic behavior 
of gn{x) for X ^ 0, 



n 

9n{x) ^ ^ 



1 



ax 



(36) 



with a = — l/[8(n- 
is reduced to 



2(n-l) 
2)(n - 1)2]. The eigen equation (31) 



a 



¥ 



C2C + ciC + Co = 



(37) 



with 

C2 = 

Cl = 

Co = 



/ 1 

a 

Vn- 1 



2na 



1 






n 1 _ 2 ^t 

4(n — 1)^ 2(n — 1) ?/^ 



+ n(/ln -£zfn) + 



n — 1 



^" 



(38) 



The eigen modes on the light line are given by the 
equation cq = 0. The exsitence of solution requires 
that Co > 0. Only the positive root with small mag- 
nitude of the above cubic polynomial will give the dis- 
persion for the modes near the light line. Since ^ is 
small, the physical solution can be approximated as 
^2 = — cq/ci — C2Co/cf — a'^CQ/cf. Those modes all have 
Ka > Xn,m for the Tn-th band. 

For n = 2, the asymptotic behavior of gn{x) is still 
given by Eq. (36), but with a non-constant coefficient 
a = [ln(x/2) + 7]/4. The eigen equation near the light 
line is reduced further from Eq. (37) to 



(/^o«)^(ln^ +7) +co 







(39) 



with Co given in Eq. (38) with n = 2. 

For the allowed eigenmodes of the 772-th hybrid band, 
one has the range < hioa < ^ri,m with (_ri,m the solutions 
of Eq. (32). The solutions near both ends of the above 
range can be obtained analytically as we have done. For 
arbitrary ^ within this range, the solution must be ob- 
tained numerically. However only when £t < —1, one 
can have eigenmodes with < Ka < Xn,i- For the m- 
th band, one has Xn^m < Ka < (,n,'m/v- Otherwise, the 
solution for the first band requires Xn,i < Ka < ^n,i/^ 
with ^^^1 > rjXn^i- The band structure for hybrid modes 
on a waveguide with e^ = — 3 and e^ = 2 is shown in Fig. 
8. The effective index of the waveguide Tieff = P/k^ is 
also evaluated and plotted in Fig. 9. 




FIG. 8: Band structure of the guided hybrid modes with n — 1 
on an anisotropic waveguide with st = — 3 and e^ = 2. Open 
circles denote the degeneracy of forward- wave and backward- 
wave modes. 
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FIG. 9: Effective index rieff as a function of the reduced radius 
k^a for the first three bands of hybrid modes with n = 1 on an 
anisotropic waveguide with s^ = — 3 and s^ = 2. Open circles 
denote the degeneracy of forward-wave and backward-wave 
modes. 



2. Energy flow of hybrid modes 

The energy flow can also be evaluated for hybrid 
modes. The expression for Sz is much more complex than 
that of the TM modes. However the final expression for 
Pz is much simpler than expected, once the integrals are 
all carried out. One has 



3 \ e'^ 1 

4^0 StKa Hca 



+ 



2n2 



1 + ^t 



2et 



{hca)^ I {Koa)'^ {tza)'^ i [' 
/3 f 1 



4fco 



Koa 



-{2gn + hn -ezfn)g'n 



2n2 



l+St 



(kqi)^ '-{Koa)^ (na)^ 



(40) 



Here /^, ^^, and h'^ are the derivatives of Z^, ^n, and /i^^ 
respectively. In this derivation, we assume et and Sz are 
real, thus the arguments of the Bessel functions are all 
real. We set A = \/gn + hn/{koaJn)- Use has also been 
made of the following integrals which can not be found 
in any mathematics manual 

Jo X Z 

rx 2 3 

/ {I'^^(x) + -Il{x))xdx = --ll{x)hUx), 

Jo X z 



/ {KHx) + -:^Kl{x))xdx = --Kl{x)g'„{xlAl) 

J ^ X z 

Explicitly, one has 






-Jn[X) + - 
X X' 



1 

X 



-xh^{x) hn{x) + ^ + -, 

2/ N 2 , , n2 1 
X9n{x) - -g„{x) - ^ - -■ 



(42) 



We point out that the above expressions for P^ can 
be readily modified for dielectric or metallic cylindrical 
waveguide with the exchange of fni'ix) = —hn{x) and 
ixf^{ix) = -xh'^{x). 

The normalized energy fiow on a waveguide with £t = 
—3 and Sz = 2 for the hybrid modes with n = 1,2,3 is 
plotted in Fig. 10. 



3. Backward- wave and forward-wave hybrid modes 

The hybrid modes have similar features as the TM 
modes that both forward waves and backward waves can 
co-exist within the same band. For koa ^ 0^ k, c^ hio^ the 
hybrid modes have d{ko(i) / d{noCi) < and are all back- 
ward waves. Only the modes near the light line can be 
forward waves. Making use of the expressions for Pz in 
Eq. (40), one can prove that dikoo) / d^noCi) > leads 
P^ > and vice versa. The degeneracy of forward- and 
backward- wave modes is located at d{koa) / d^Koa) = or 

(Pz) = 0. 

For the hybrid modes with n = 1,2, the modes near 
the light line are forward waves since dikoo) / d{no(i) > 0. 
However as can be seen in Fig. 10, the whole first band 
of the hybrid modes with n = 3 are backward waves. 
Hybrid modes with higher n will have more bands to be 
all backward waves. If one further increases the angular 
index n of the hybrid modes, more hybrid mode bands 
will be all backward waves. 
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FIG. 10: (Color online) Normalized energy flow (Pz) for the 
first few bands of the hybrid modes with n = 1,2,3 on an 
anisotropic waveguide with £t = — 3 and £2 = 2. Open circles 
denote the degeneracy points where (Pz) = 0. 



III. SLOW AND TRAPPED LIGHT 

Recently Tsakmakidis et at [31] proposed to trap light 
in a tapered waveguide with double negative index. The 
indefinite index waveguide we have studied so far in this 
paper can also be used to slow down and trap light. 
These waveguides can thus be used as optical buffers [32] . 
The reason is that unlike the ordinary optical fiber, these 
waveguides support both forward and backward waves. 

For the anisotropic waveguide we have considered, 
St < 0, one has P]" < and P^°^* > if one sets /5 > 0. If 
Pz = P^^-\-P^^^ < 0, the mode is a backward mode since 
the total energy fiow is opposite to the phase velocity. 
Otherwise, the mode is a forward mode. At the criti- 
cal radius ac^ the backward and forward modes become 
degenerate, the energy fiow inside the waveguide cancels 
out that in the air. One can prove that at the critical ra- 
dius ac where P^ = 0, the group velocity is indeed zero. 
One does not need to know the material dispersion to 
locate the zero group velocity point. This is due to the 
fact that for these waveguides, the dispersion due to geo- 
metric confinement dominates the material dispersion at 
and around the critical radius. 

The unique properties of the modes on anisotropic 
waveguide can be used to slow down and even trap light. 
Even though the waveguide supports infinite number of 
both TM and hybrid modes at any fixed radius and fre- 
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quency, with appropriate laser coupling, the excitation 
of the hybrid modes in the waveguide can be suppressed 
or even eliminated. Among the TM modes, the first TM 
mode will be more favorably excited. Furthermore, due 
to the material dissipation, the first TM mode will propa- 
gate the longest distance. The rest of the TM modes will 
all decay out at about half the decay length of the first 
TM mode. It is the first TM band which can be used for 
slow light application. Unlike the double negative waveg- 
uide [31], the anisotropic waveguide will slow down and 
trap light if one increases the radius to the critical radius. 
A sketch of a slow light waveguide is shown in Fig. 11. 




FIG. 11: (Color online) A sketch of the slow light waveguide. 



IV. REALIZATION OF EXTREMELY 
ANISOTROPIC NANOWIRES 

These extremely anisotropic media can be realized in 
metamaterials. According to the effective medium theory 
[26, 33], one has for multilayered structure of dielectric 
Sa and metal e^n, the effective permittivities are 



St 



fSm + (1 - f)ea, 
fSa + (1 - /)£m ' 



(43) 



Here / is the filling ratio of the metal. For / > /min = 
^a/i^a — 3?^m), ouc has ^Sf < 0. A realization of the 
anisotropic nanowire is shown in Fig. 12. 



lattice spacing D 




FIG. 12: (Color online) A sketch of the realization of an 
anisotropic nanowire made of alternative disks of metal and 
dielectric. 



We first consider a metamaterial waveguide at a fixed 
wavelength. For silver at A = 488 nm, one has Em = 
— 9.121 + 0.304z [34]. A nanowire made of aler native disks 
of silver and glass {sa = 2.25) of equal thickness will have 
£t = -3.436 + 0.152zand6^ = 5.971 + 0.065i by using the 
effective medium theory. Here the disk thickness is 10 nm 
for both materials. For example if one sets a = 60 nm, 
one has koa = 0.7725. The first three TM modes will 
have (3a = 2.2525 - 0.0747^,4.9318 - 0.1419^,7.4031 - 
0.2088i. Thus one has Xp ~ 167 nm and phase refrac- 
tive index neff = 2.92 for the first TM mode. The decay 
length is 803, 423, and 287 nm, respectively. After trav- 
eling about 420 nm along the nanowire, only the first one 
will survive. 

Finite-difference time-domain (FDTD) simulations 
[35] were performed to obtain the effective index neff of 
the metamaterial nanowire. The procedure is the follow- 
ing. We illuminate the free-standing nanowire of finite 
length with a Gaussian beam, then get Ez after the ter- 
mination of the simulation. The length of the waveg- 
uide is set to be larger than the decay length of the first 
TM mode. We get the phase from Ez, then determine 
p. Though the waveguide supports infinite number of 
modes including TM and hybrid modes, our method is 
legitimate due to the following two reasons. First that 
the excitation of hybrid modes is small due to the pro- 
file of the incident Gaussion beam. So mainly the TM 
modes are excited. Second that due to the dissipation in 
the metal, after certain distance, only the first TM mode 
will survive. Thus the phase propagation is mainly due 
to the first TM mode. The amplitide and phase propa- 
gation of Ez along the above metamaterial nanowire is 
shown in Fig. 13. 

Abs(E^) 
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FIG. 13: (Color online) FDTD simulation of the amplitude 
and phase propapation of the longitudinal electric field Ez 
along the nanowire with radius a = 60 nm at A = 488 nm. 
The metamaterial nanowire consists of alternative disks of 
silver and glass disks of thickness 10 nm. 

The relation between the effective index neff and the 
nanowire radius a is shown in Fig. 14. Very good agree- 
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merit between FDTD simulations and analytical results 
has been obtained. However for small radius, there is 
noticeable discrepancy. This is expected since when the 
radius is comparable with the lattice spacing of the mul- 
tilayered metamaterial, the effective medium theory will 
fail. We have also performed FDTD simulations for the 
nanowire with smaller lattice spacing. Better agreement 
is indeed obtained. 
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FIG. 14: (Color online) The effective index rieff of the first 
TM band on a nanowire with different radius at A = 488 nm. 
The nanowire is made of alternative disks of silver and (see 
Fig. 12). The disk thickness is 10 nm. Filled circle (green) is 
obtained from FDTD simulations. The dashed line (blue) is 
the fitting of simulation data. The solid line (red) is calculated 
from band equation with effective index £t — —3.436 + 0.1 52i 
and e^ = 5.971 + 0.065i. 

We also consider the band structure for different fre- 
quencies. The permittivity given by Eq. (24) can be 
realized through the multilayered heterostructure with 
Drude metal with £^ = 1 — kp/ko and dielectric Sa- The 
nanowire is made of alternative disks of a Drude metal 
and a dielectric. The band structure and the effective in- 
dex neff of the TM modes are shown in Fig. 4 and Fig. 5, 
respectively. One noticeable feature of these bands is the 
flatness of each band, which indicates small group veloc- 
ity. We have erformed the FDTD simulation for different 
frequencies for nanowire with a fixed radius. The results 
are shown in Fig. 15. Again good agreement between 
FDTD simulations and analytical results is achieved. 



V. CONCLUSIONS 

Indefinite index materials can be used to achieve neg- 
ative refraction [25] and hyperlensing [27, 28]. They can 
also be used as superlens [26]. In this paper, we con- 
sider the wave propagation along a cylindrical waveguide 
with anisotropic optical constant. We have derived the 
eigenmodes equation and obtained the solutions for all 
the propagation modes. The field profiles and the en- 
ergy flow on the waveguide are also analyzed. Closed- 
form expressions for the energy flow for all the modes 




0.5 
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FIG. 15: (Color online) The effective index ries for the TM 
modes on a nanowire with radius a = 40 nm. The nanowire 
is a stack of equally thick alternative disks made of a Drude 
metal Sm = 1 — kp/ko(ko + iT) and glass. Here kpa — 1.64 
and Fa = 0.0155. The analytical curves (solid) are calculated 
by using real Sm- 



are derived. For extremely anisotropic cylinder where 
the transverse component of the permittivity is negative 
and the longitudinal is positive {st < 0, e^ > 0), the 
waveguide supports TM and hybrid modes but not the 
TE modes. Among the supported TM modes, at most 
only one mode can be forward wave. The rest of them 
are backward waves. 

The case that St > ^ and e^ < can be discussed simi- 
larly. Anisotropic waveguides of cross section other than 
circle were also considered. The results will be published 
elsewhere. 

Possible realization of these extremely anisotropic 
nanowires are proposed. Extensive FDTD simulations 
have been performed and confirmed our analytical re- 
sults. 

Two unique properties have been revealed for the 
modes on nanowire waveguides made of indefinite in- 
dex metamaterials. The first is that the backward- 
wave modes can have very large effective index. These 
nanowires can be used as phase shifters and filters in 
optics and telecommunication. The second is that the 
waveguide supports modes of zero group velocity. This is 
due to the fact that the waveguide can support both for- 
ward and backward waves at a fixed radius. If the waveg- 
uide is tapered, at certain critical radius, the two modes 
will be degenerate and carry zero net energy flow. At 
other radii, these waveguides support modes with small 
group velocity. These waveguides can also be used as 
ultra-compact optical buffer [32] in integrated optical cir- 
cuits. 
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One thus has 



APPENDIX A: PADE APPROXIMANT OF fn{x) 



fn,j = (4 - j){jA - 4b)/[jc{jA - 4a). We further define 

Cfl = 3/n,2//n,l, 0-2 = 3/n,3//n,2. (A3) 



Consider the function /^(x) = J!^{x)/[xJn{x)]. Let ^^^r some manipulations of the algebra, one obtains 



be the m-th zero of the Bessel function Jn{x). We 
also denote x^ ^ as the zeros of Jn{x). In order to get 
the eigen modes on the anisotropic waveguide easily, we 
may need the inverse function f^^ix) in the interval 
[xn.m^x'^rn]' We cousidcr the Fade approximant to the 
function fn{x)^ 



fn{x) 



{X 



^){^ 



h) 



c{x ■ 



^){X■ 



a) 



(Al) 



Instead of fixing the three unknowns through the coef- 
ficients of the Taylor expansion of fn{x)^ here we deter- 
mine them by the exact values of fn{x) at some points. 
Since there are three unknowns, we only need the value of 
fn{x) at three points. For simplicity, we evaluate fn{x) 
at three evenly spaced points 



Jn 



Jn\Xri 



-iA/4) 



(A2) 



4V cri+(Trl-2/' 



(Tl -r ^2 

1 A -26 

'hA-2a 



(A4) 



The inverse of the function fn{x) can be obtained as one 
of the roots of a quadratic equation. As it turns out, the 
expression obtained in this way gives very good approx- 
imation to f~^{x). 

Once the inverse function is obtained, one has Ka = 
fo^{y) with y = go{i^oCi)/£z foi" the TM modes. Together 
with the asymptotic expressions of ^n(^) and hn{x), Eq. 
(Al) can also be used to obtain approximate solutions of 
the hybrid modes. 
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